Multiprocessor Approximate Message Passing with Column-Wise Partitioning

Solving a large-scale regularized linear inverse problem using multiple processors is important in various real-world applications due to the limitations of individual processors and constraints on data sharing policies. This paper focuses on the setting where the matrix is partitioned column-wise. We extend the algorithmic framework and the theoretical analysis of approximate message passing (AMP), an iterative algorithm for solving linear inverse problems, whose asymptotic dynamics are characterized by state evolution (SE). In particular, we show that column-wise multiprocessor AMP (C-MP-AMP) obeys an SE under the same assumptions when the SE for AMP holds. The SE results imply that (i) the SE of C-MP-AMP converges to a state that is no worse than that of AMP and (ii) the asymptotic dynamics of C-MP-AMP and AMP can be identical. Moreover, for a setting that is not covered by SE, numerical results show that damping can improve the convergence performance of C-MP-AMP.Comment: This document contains complete details of the previous version (i.e., arXiv:1701.02578v1), which was accepted for publication in ICASSP 201

Pulsar Timing Probes of Primordial Black Holes and Subhalos

Pulsars act as accurate clocks, sensitive to gravitational redshift and acceleration induced by transiting clumps of matter. We study the sensitivity of pulsar timing arrays (PTAs) to single transiting compact objects, focusing on primordial black holes and compact subhalos in the mass range from $10^{-12} M _{\odot}$ to well above $100~M_\odot$. We find that the Square Kilometer Array can constrain such objects to be a subdominant component of the dark matter over this entire mass range, with sensitivity to a dark matter sub-component reaching the sub-percent level over significant parts of this range. We also find that PTAs offer an opportunity to probe substantially less dense objects than lensing because of the large effective radius over which such objects can be observed, and we quantify the subhalo concentration parameters which can be constrained.Comment: 18 pages, 6 figure

Shedding light on a living lab: the CLEF NEWSREEL open recommendation platform

In the CLEF NEWSREEL lab, participants are invited to evaluate news recommendation techniques in real-time by providing news recommendations to actual users that visit commercial news portals to satisfy their information needs. A central role within this lab is the communication between participants and the users. This is enabled by The Open Recommendation Platform (ORP), a web-based platform which distributes users' impressions of news articles to the participants and returns their recommendations to the readers. In this demo, we illustrate the platform and show how requests are handled to provide relevant news articles in real-time

Transformations of Boolean Functions

Boolean functions are characterized by the unique structure of their solution space. Some properties of the solution space, such as the possible existence of a solution, are well sought after but difficult to obtain. To better reason about such properties, we define transformations as functions that change one Boolean function to another while maintaining some properties of the solution space. We explore transformations of Boolean functions, compactly described as Boolean formulas, where the property is to maintain is the number of solutions in the solution spaces. We first discuss general characteristics of such transformations. Next, we reason about the computational complexity of transforming one Boolean formula to another. Finally, we demonstrate the versatility of transformations by extensively discussing transformations of Boolean formulas to "blocks," which are solution spaces in which the set of solutions makes a prefix of the solution space under a lexicographic order of the variables

Markov chain Monte Carlo and expectation maximization approaches for estimation of haplotype frequencies for multiply infected human blood samples

Background Haplotypes are important in anti-malarial drug resistance because genes encoding drug resistance may accumulate mutations at several codons in the same gene, each mutation increasing the level of drug resistance and, possibly, reducing the metabolic costs of previous mutation. Patients often have two or more haplotypes in their blood sample which may make it impossible to identify exactly which haplotypes they carry, and hence to measure the type and frequency of resistant haplotypes in the malaria population. Results This study presents two novel statistical methods expectation–maximization (EM) and Markov chain Monte Carlo (MCMC) algorithms to investigate this issue. The performance of the algorithms is evaluated on simulated datasets consisting of patient blood characterized by their multiplicity of infection (MOI) and malaria genotype. The datasets are generated using different resistance allele frequencies (RAF) at each single nucleotide polymorphisms (SNPs) and different limit of detection (LoD) of the SNPs and the MOI. The EM and the MCMC algorithm are validated and appear more accurate, faster and slightly less affected by LoD of the SNPs and the MOI compared to previous related statistical approaches. Conclusions The EM and the MCMC algorithms perform well when analysing malaria genetic data obtained from infected human blood samples. The results are robust to genotyping errors caused by LoDs and function well even in the absence of MOI data on individual patients

Twist instability in strongly correlated carbon nanotubes

We show that strong Luttinger correlations of the electron liquid in armchair carbon nanotubes lead to a significant enhancement of the onset temperature of the putative twist Peierls instability. The instability results in a spontaneous uniform twist deformation of the lattice at low temperatures, and a gapped ground state. Depending on values of the coupling constants the umklapp electron scattering processes can assist or compete with the twist instability. In case of the competition the umklapp processes win in wide tubes. In narrow tubes the outcome of the competition depends on the relative strength of the e-e and e-ph backscattering. Our estimates show that the twist instability may be realized in free standing (5,5) tubes.Comment: 4 pages, 1 figur

Spectroscopic Constants and Line Positions for TiO Singlet States

consistent set of spectroscopic constants for the a1Δ,d1Σ+,b1Π,c1Φ, and f1Δ states of 48Ti16O has been determined from analysis of the b1Π–a1Δ,b1Π–d1Σ+,c1Φ–a1Δ, and f1Δ–a1Δ systems. Three Fourier transform emission spectra have been used for the analysis. New bands of the b1Π–a1Δ and c1Φ–a1Δ systems have been fitted. The first analysis of the c1Φ–a1Δ system using Fourier transform spectra is also provided. Extensive and improved line positions are measured. TiO is prominent in the spectra of oxygen-rich cool stellar objects and may be present in hot-Jupiter exoplanet atmospheres

Line Lists for LiF and LiCl in the X1Σ+ Ground State

Vibration–rotation line lists for 6LiF, 7LiF, 6Li35Cl, 6Li37Cl, 7Li35Cl, and 7Li37Cl in the X1Σ+ ground states have been prepared. The rovibrational energy levels have been calculated using potential energy surfaces determined by direct potential-fitting employing the rotational and rovibrational transition frequencies of all isotopologues, and required the inclusion of Born–Oppenheimer breakdown terms. Dipole moment functions calculated ab initio at the MRCI/aug-cc-pwCV5Z level have been used for line strength calculations. Partition functions for temperatures up to 5000 K have been calculated. LiF and LiCl are predicted to be present in the atmospheres of hot rocky exoplanets, brown dwarfs, and cool stars

Cosmological Tension of Ultralight Axion Dark Matter and its Solutions

A number of proposed and ongoing experiments search for axion dark matter with a mass nearing the limit set by small scale structure (${\cal O} ( 10 ^{ - 21 } {\rm eV} )$). We consider the late universe cosmology of these models, showing that requiring the axion to have a matter-power spectrum that matches that of cold dark matter constrains the magnitude of the axion couplings to the visible sector. Comparing these limits to current and future experimental efforts, we find that many searches require axions with an abnormally large coupling to Standard Model fields, independently of how the axion was populated in the early universe. We survey mechanisms that can alleviate the bounds, namely, the introduction of large charges, various forms of kinetic mixing, a clockwork structure, and imposing a discrete symmetry. We provide an explicit model for each case and explore their phenomenology and viability to produce detectable ultralight axion dark matter.Comment: 11 pages, 3 figure